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Abstract. We study the equivalence of two — order-by-order Einstein’s equation and Re¬ 
duced action — approaches to cosmological perturbation theory at all orders for different 
models of inflation. We point out a crucial consistency check which we refer to as ‘Con¬ 
straint consistency’ condition that needs to be satisfied in order for the two approaches to 
lead to identical single variable equation of motion. The method we propose here is quick 
and efficient to check the consistency for any model including modified gravity models. Our 
analysis points out an important feature which is crucial for inflationary model building i.e., 
all ‘constraint’ inconsistent models have higher order Ostrogradsky’s instabilities but the re¬ 
verse is not true. In other words, one can have models with constraint Lapse function and 
Shift vector, though it may have Ostrogradsky’s instabilities. We also obtain single variable 
equation for non-canonical scalar held in the limit of power-law inhation for the second-order 
perturbed variables. 


1 Introduction 


Inflation has now become an integral part of the standard model, that can eliminate cosmolog¬ 
ical initial value problems, explain homogeneities as well as inhomogeneities and observation 
of anisotropic Cosmic Microwave Background Radiation (CMBR)[1]. It is a period of accel¬ 
erated expansion in the very early universe and it occurs around 10^^ GeV which is much 
remote in time compared to the terrestrial experiments. Inflationary cosmology has two key 
theoretical aspects. One is the approximation schemes employed in solving gravity equations. 
The other is the inflationary model building inspired by particle physics or a fundamental 
theory of Quantum gravity. The problem with any theory of gravity is that it is typically 
highly non-linear, so one has to rely on approximation schemes to match the observations. 
Primarily there exist two formalisms to deal with the non-linear equations: 

• The separate universe approximations [2-4] with either gradient expansion theory or 
AN formalism. [5-11] 

• Gauge invariant cosmological perturbation theory[12-34]. 

The temperature fluctuations as observed in CMB is ~ 10“®, hence it is consistent to use 
order-by-order perturbation theory to match with observations[12, 35, 36]. In the first order, 
one assumes that the perturbed fields are linear. This implies that the 3-point and higher or¬ 
der correlation functions are zero. In the second order, the interactions of the first order need 
to be included, hence, leading to non-zero 3-point functions. Also it is widely believed that 
the detection of these 3-point correlation functions can reduce the field space of inflationary 
models[14, 17, 37]. 

With respect to inflationary model building, the proposed theories are primarily pre¬ 
ferred through simplicity. In the case of the canonical scalar held, the simplest, 60 e-foldings 
of inflation require the potential to be flat, which is in contradiction with particle physics 
models[38, 39]. Non-canonical scalar field model[40-42] removes the dependence of the poten¬ 
tial, however it leads to time dependence of the speed of perturbations and makes it difficult 
to be compared with CMB observations [43]. In order to seek more generalized fields, scalar 
fields with higher time derivatives in action are considered[44-47]. Beside these, modified 
gravity models, specifically f{R) lead to accelerated expansion in the early universe. 

There are two mathematical procedures that are currently used in the literature to 
study gauge invariant cosmological perturbation theory: Hamiltonian formulation (ADM 
formulation) [34] and Lagrangian formulation[12-33]. Since gravity and matter are coupled 
to each other, one can write the full action and vary the action with respect to metric and 
matter fields to obtain general equations of motion (e.g., Einstein’s equation in General rel¬ 
ativity). Those equations can be expanded in terms of perturbed variables (metric and field 
variables) and one can write down equations in the perturbation theory [48]. In the action 
formalism [49], the action is expanded to the required perturbed order in terms of the per¬ 
turbed variables (metric and field variables) and then varied the action with respect to these 
variables. For example, to obtain perturbed equations in the first order, we need to expand 
the action to second order and vary the action with respect to the first order perturbed vari¬ 
ables. This formalism can be extended to obtain perturbed equations of motion up to any 
order. In the reduced action formalism, constraint variables are replaced in the action using 
constraint equations so that we can rewrite the action only in terms of dynamical variables. 

Since the matter fields (Non-canonical & Galilean scalar field model) and Gravity are 
highly non-linear, it is not clear whether the two approaches, i.e., Einstein’s equations writing 
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in order-by-order perturbation theory and action/reduced action formalism, lead to the same 
equations of motion. In Ref. [50], it was shown that when the metric perturbations are frozen 
then the two approaches do not, in general, lead to the same expressions. In this work we 
address the issue by including the metric perturbations in the theory. 

In the next section, we study higher order cosmological perturbation theory for a single 
scalar field minimally coupled to gravity and show the equivalence of the two approaches at 
all orders. We point out a crucial and novel consistency check which we refer to as ‘constraint 
consistency’ condition that needs to be verified. We also show that this provides a fast and 
efficient way to check the consistency and apply it to minimally coupled non-canonical scalar 
field. 

In section 3, we apply the ‘constraint consistency’ condition to many inflationary mod¬ 
els that are proposed in the literature. First we check the theory with higher derivative 
Lagrangian models minimally coupled to gravity. Then we extend the procedure to other 
different types of models like modified gravity models and modihed gravity with higher order 
matter Lagrangian. Appendix A contains some of the derived expressions used in section 2 
and in Appendix B, we obtain a single variable equation of motion for non-canonical scalar 
fields in terms of second order perturbed variables. 

In this work, the number of space-time dimensions is 4 and the metric signature we use 
is [—, , K = SttG, c = 1. 

2 Consistency of Higher order pertnrbations in two different approaches 

The action for gravity sourced by a single, non-minimally coupled scalar field {(p) is, 

S = j + ( 2 . 1 ) 


where 

£^ = PiX,p) + G{X,p)Dp, X = ^g^''d^pd^ip, □ =(2.2) 

is the Lagrangian for the Galilean field which is the most general scalar field model leading to 
second order equations of motion. Varying the action with respect to metric gives Einstein’s 
equation, 

R/iu - R = K (2.3) 

where the stress tensor is, 

= gf,u{P + Gxg'^^daXdpip + G^g^'^dapdpp} ^ 

- {Px + 2G^ + Gxapjd^ifdyif - 2Gxdf,Xdyp 

For simplicity and to obtain the physical features, we consider only single scalar field 
theory minimally coupled to the gravity. In Section 3.4, we look at modified gravity models. 
Variation of the action (2.1) with respect to the scalar field leads to the following equation 
of motion, 

{2G^ — 2XGx<p + Px}G\ip — {Pxx + 2Gx^}d^pd^ip — 2X{Gip -|- Px<p} + Pip 
- GxW,pup':'' - {UpY + R^yd^^ipd’^ip} - Gxx{d^,xd>^x + {d^pd>^x}ap}} = o 
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As one can see, although the Lagrangian is of the form, Cm = Cm{c^ d‘^tp, t), i.e., it 
contains higher time derivatives of the scalar field but equations of motion are second order, 
thus does not suffer from Ostrogradsky’s instability [51]. With G{X, ip) = 0 the field becomes 
non-canonical. Further fixing P = —X — V{(p), where V((p) is the potential, the Lagrangian 
corresponds to canonical scalar field. 

The four-dimensional line element in the ADM form is given by, 

= g^udx^dx'^ 

= -{N^ - NiN^)djf + 2Nidx^dr] + jijdx^dx^ ( 2 . 6 ) 


where N(x^) and Ni{x^) are Lapse function and Shift vector respectively, 'jij is the 3-D space 
metric. Note that, in the case of Galilean model, N(x^) and Ni(x^) are the gauge constraints 
and variation of action (2.1) with respect to those lead to Hamiltonian and Momentum 
constraints, respectively. 

Action (2.1) for the line element (2.6) takes the form, 

S = J d^xV^{^ ( 2 . 7 ) 


where Kij is extrinsic curvature tensor and is given by 


- 2N 
K = 


do^ij ^i\j ^j\i 




Perturbatively expanding the metric and the scalar field about the flat FRW spacetime, 
we get. 


<700 = —0^(1 + 2e(/>i -|- 4>2 + •••) (2-8) 

9oi = Nj = a‘^{€diBi + -e^diB2 + ■■■) (2.9) 

9ij = ~ 2e^i — e^t />2 — + 2eEmj + E2\ij + ...} (2-10) 

ip = ipo + epi +-e^(p2 +■■■ (2.11) 

where ‘e’ denotes the order of the perturbation. Note that we have ignored the vector and 
tensor part of the metric perturbations. Although in the first order, the scalar, vector and 
tensor perturbations decouple, the three types of perturbations are coupled in higher order. 
We assume that the vector and tensor contributions are small and can be neglected at all 
orders. 

To determine the dynamics at every order, we need five scalar functions {(f), B, if), E 
and p) at each order. Since there are two gauge choices, one can fix two of the five scalar 
functions. In this work, we choose flat-slicing gauge, i.e., 'i/i = 0, Fi = 0 at all orders. 


9oo — —a^(l + 2e(/ii -|- (f)2 + ...) 

(2.12) 

9oi ^ Ni = a^iediBi + ^e^diB2 + ...) 

(2.13) 

9^J = 

(2.14) 

1 2 

if = {p^^eipi + -e ip 2 + ... 

(2.15) 
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In the next subsection, we obtain the equation of motion of the second order perturbed 
quantity in single variable form for non-canonical scalar field using order-by-order perturbed 
Einstein’s equation. In subsection 2.2, we use reduced action approach. To confirm or 
infirm the result of Ref. [50], that the two approaches lead to different results we focus on 
non-canonical scalar field, i.e., setting G(X,ip) = 0. 

2.1 Order-by-Order Einstein’s equation approach 

For the background, = diag{—a?‘^a?‘,a?,a?‘), equations (2.3) and (2.5) lead to, 

— {Px^'o"+ Pa^)=n^ (2.16) 

-2—+n^ = KPa^ (2.17) 

a 

PxVo - Pxx^Q^Q + Pxip'^Q + 2Px<7’o'R + PxxP^'q = 0 (2.18) 

Equations (2.16) and (2.17) are zeroth order 0-0 and i-j Einstein’s equations where as 
equation (2.18) is the zeroth order equation of motion of the scalar field. Similarly, the first 
order 0-0, 0-i Einstein’s equations and equation of motion of the perturbed scalar field are, 

HV^Bi = -{Px4>i‘^'q + 2Pa^(j)i -|- + Pxx4>i^o‘^0'~‘^ - Pxxy^Wo^a~‘^ + 

PxipV>'o^^i + P^^ia^) (2.19) 

P(l>i = ( 2 . 20 ) 


- Px<^ia^ - Pxx4>i‘fo^ + Pxx^Wo^ - Pxx^4>Wo‘^ + PxXif^W^ “ P>p(f>ia‘^ 

- P^^a^^pi + Px(t>iVQa^ + - 2Px(p[Pa^ - 

+ ^Pxx‘pWi‘fo + PxXif’fWo^'fi ~ Pxip^Wi^'^ ~ Pxip>^oa^‘Pi - Px<fi<pv'(i^ 

+ 2Px<?^i7^o^a^ + Px^'o^^Bio^ + Pxx(t>iP^'^ — Pxx^'iH^o^ — Pxxx(t>iP^'^ oT"^ 
+ Pxxx(t>Wov'o^o.~‘^ + Pxxxv'iH'^'o^a~‘^ — Pxxx^Wo^o^o,~'^ — Pxx^P<f'{^<Pi 

- 2Px!p^QP‘Pia^ = 0 


( 2 . 21 ) 


Note that, there are no </>'/ and B” terms in the above three equations and equations 
(2.19) and (2.20) are, as expected, the constraint equations corresponding to Lapse function 
and Shift vector. Hence, ())i and Bi are constraints and we can eliminate them from the first 
order equation of motion of the scalar field (2.21). In first order, single variable equation for 
non-canonical scalar field in terms of Mukhanov-Sasaki variable (v) is. 


where. 


II 2x72 ^ n 

V — c^V V - V = 0 


_ ^ 2 ^ Px 

P ’ ^ Px + 2XPxx 


( 2 . 22 ) 


(2.23) 
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Similarly, perturbed second order 0-0 and 0-i Einstein’s equations for non-canonical 
scalar fields at second order are, 

- - 2S^md,Bidj(l)i - -t S’^^d.kBidjiBi 

+ k{-^Px 5"^diBidjBitp'o^ - PS^^diBidjBio^ - ]^Px(t>2^'o^ - P(t>2a^ + 2Px4>‘iipQ^ + iPfj^la^ 

7 4 2 1 

-I- 2Pjc<^i<Po‘/?i + i^Pxx(t>Wo a~‘^ - ^>Pxx4’ivWo + Pxv4‘i‘P'o Vi - 2^x‘fW2 

— ^Px^'i^ — 2Pxx4‘2y^o'^a~‘^ — -PxxS’'^ diBidjBi(f'Q^a~‘^ + 2Pxx^'o^^'i^a~'^ + ^Pxx^Wf^ 

- PxxS"^diBidjipnpQ^a~‘^ - ]^Pxx5"^a~‘^ - Px^^’o^pWi - ]^PxvV>'o^V>-2 

— 2^xxx(t>\^'{)^oT'^ + Pxxx4’wWo^~ ^Pxxx^'o^'^'^a~‘^ — ^Pxvv^p'o^^fi — Pxxip4>i^o‘^a~^(pi 
+ PxXipfWo^(^~^^i - ^PxS'^diifidjipi - ^P^(p2a'^ - ^P^v^pW) = 0 


(2.24) 


- 25^^ndjBid,kBi + -t 4.(t>iHd^(j)i - Hd^(t)2 - 6^’^dj(t)id^kBi 


+ K{--PxPo^iP2 - Pxp'idiipi - Pxxfl^idiipiPo a ^ + Pxxp'id.ipiip^ a ^ - PxipP'od.ipiipi) = 0 

(2.25) 

and the equation of motion of the scalar field is, 


CxPx + CxxPxX + CxXxPxXX + CXXXXPXXXX + CxXXifiPxXXif 
+ CxxifPxxip + CxxififiPxx^pip + CxipPxifi + CxipipPxifiip + CxipipipPxipif^ (2.26) 

4“ CipPi^ 4“ C(fipPifip 4“ Ci^ipipPi^ipi^ — 0 

where Cx-, Cxx, ■■■ are all second order perturbed quantities and Px,Pxx, Pxip, ■■■ are back¬ 
ground quantities. The explicit form of C’s are given in Appendix A. 

It is important to note that the second order equations also do not contain (p'^ and/or 
^2 • Hence one can obtain a single variable equation of motion of non-canonical scalar field 
at second order. Malik et al obtained the single variable equation of motion in second order 
for canonical scalar field[52]. Also note that V’ evaluated in the flat-slicing gauge is a gauge 
invariant quantity and directly related to comoving curvature perturbation 7^/curvature per¬ 
turbation on uniform-density hypersurfaces (/ [53]. 

2.2 Reduced Action approach 

In the reduced action approach, which is now a popular way to calculate non-gaussianity, 
one perturbs the field variables (p) in the action and expands the action to the required 
order. In other words, one assumes a priori the form of the metric and the matter variables in 
the lowest order and expands order-by-order. For instance, in the case of FRW background, 
the action (2.7) becomes. 


(0)5^^ = I d\i Pa^ - 3-a 

K 


1 /2 


(2.27) 


Varying the above action with respect to metric variable a{r]) and p^oiv) leads to the 
equations (2.17) and (2.18). Note that, as expected, these two equations are independent of 
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each other since a{rj) and (fioiv) dynamical variables. To obtain first order (in e) equations, 
one expands action (2.7) upto second order of e. In general, varying the order action with 
respect to order perturbed variables leads to (n — order perturbed equations. It may 
be worth noting that a given order equations of motion can be obtained in several ways, e.g., 
varying first order action with respect to first order variables leads to zeroth order equations 
of motion. 

Expanding the action (2.7) to the second order, only in terms of first order variables 
we get 


i2)gNC ^ J - PxfWo'a^ + Pxfiy^WW - \Pxg^'i 


2 2 
a - 


Px6"^(Pod,Bidjipia‘^ + ^PxS’-^diipidjipia'^ + ^Pxxfly^’o^ - Pxxfiy^Wo^ + ^Pxx^'o^ 
^Py,^(pW + Px^fiv'o^a^pi - Pxip'p'o'pW'Pi + Py,(t>ia"^<pi + ^P5"^d^BidjBia'^ - ^Pcjfia^- 

AC ^ AC ^ AC ^ AC 

Z AC 

(2.28) 


After integrating by-parts, and dropping off boundary terms, we get. 


{2)gNC ^ 


J d"^x(^ 


1 „ „ „ „ ,2 2 n ^ 2 , / 2„2 , n i ^ 


-Pxd"^^d^BidjBifi^'"of - Pxfl^'oa^ P PxfiPoPW “ d^x'p['^a^+ 


PxS^'^PodtBidjpia^ + ^PxS'’^diipidjpio^ + ^Pxxflp'o^ - PxxfipWo^ + ^Pxxp'o^fii^+ 


+ Px^fip'o^a'^pi - Pxipp'oPWpi + Py.fia'^Pi + ^P6^^d.BidjBia'^- 

IpfW - ^fiS^^-a'd.^Bia + -d.Bid^Bia'^ - ^ 

Z AC ^ AC ^ AC ' 

(2.29) 


Varying action (2.29) with respect to ipi, we obtain first order equation of motion of the 
scalar field same as (2.21). Similarly, varying action with respect to fii and Bi gives same 
equations as (2.19) and (2.20) respectively, i.e.. 


^1 

6(pi } 


6S2^ 
64>1 / 


SBi J ipi,4>i 


= 1^* order Equation of motion of the scalar field 
= 1®* order Hamiltonian constraint 
= 1®* order Momentum constraint 


Similarly, we can expand (2.7) upto fourth order by expanding the field variables 
^ 2 , B 2 ) and vary the action with respect to second order perturbed field variables to 
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obtain second order equations. Fourth order action containing only ip 2 terms are, 


= j d‘^x{^Px(l>2‘fWW + ^Px&''WW2diBidjBia^ - '^Px‘fW2(l>W + - ^Bx‘P2^a^ 

-Pxd'‘^‘PodiB2djip2a^ — -Px<l>iS'^ ‘fiodiBidj(p2a,^ + -Px^''^ v'idiB\dj^p20,^ + -PxS''^ ‘f2diBidjipia?+ 

-Px&P diip2djip20? + -Pxx‘P24‘i‘fo^ — ‘^Pxx<l>iv'i'P2‘P'o^ + -Px X diBidjip2‘Po^ + -Pxx4>i^'‘^ diLp\djip2ip'Q^ 

o 2 11 

-^Pxx4>2<P2‘p'o^ — -PxxSP(p'2diBidjBnpQ^ + -Pxx‘Po‘P2Vi^ - -Pxx^’'^‘PidiBidjip2‘Po^- 

-Pxx^''^ v'ov'idnpidjtp2 + -Pxx‘p'o^V2^ — -PxxS''^ V2diBidjipnpQ^ — -PxxS''^ Vo'P2di‘Pidj‘Pi + 

1 o 1 4 

- P,p,fiifi20.‘^ + -Pxip<p2‘fo^a'^‘P2 + -Px<fiS'‘^ diBT_djB\ipQ^a'^ip2 — -Px,p4‘i'p'o^o?V2 + -Pxip4>i'Po‘P'i’^^‘P2 + 
o 4 4 11 

-Pxip<t>iVo‘P2a'^ipi — -Px^‘fio‘fi20''^‘f2 — -Px^‘p'i^a'^‘f2 — -Px<pV’i¥’2“^¥’i + -PxipS^^ ‘fodiBidjipia'^(fi2+ 
-PxipS^^‘fodiBidjtf2a^<fii + -Pxip&^^diipidjipia^if2 + -Pxipi^^dnp\djip2a,^‘fi - -Pxxx<fi24'i‘fQ^O'~'^+ 
-PxXX<I>1‘pWWq^°’~'^ - -PxXX<P 2 ‘Pq^'-Pi^<^~^ + -Piptp<p‘fl 0 .'^'fi 2 + -PxXip4'l‘fo*‘P2 - -Pxx,p<l>iV>Wo^‘P 2 - 
-PxX,p4>lV>Wo^Vl + ■^PxX^'f'o^‘f'l^'-P2 + -PxX^vWWo^'fi + -Px^ip<Pl‘fo‘^a.'^‘Pl‘f2 - -Pxip<p‘pWl°^‘fi‘P2- 


^Pxip,p‘pW2vio-'^ + ^PxfpiS'^ di(fidjifi2a'^ + ]^Pipip4>ia^‘Pi‘P2 + ^Pip<t>2a^V2 + ^P^pS^^ diBidjBia^ifi2 
- jPp<l>W‘P2) 


(2.30) 


Varying the action with respect to (p 2 leads to the same equation of motion of (p 2 
(2.26). Similarly second order equations of (f >2 and B 2 can be obtained from varying fourth 
order action with respect to (/)2 and B 2 (2.24) and (2.25), respectively. This mechanism can 
be extended up to any order and we can generalize that, equations obtained from both the 
approaches are identical and there are no ambiguities as discussed in Ref. [50]. 

Another way of seeing constraints is, action (2.29) or (2.30) contain no time derivative 
of (/)i, i?i, </>2 and B 2 , i.e.. Lapse function and Shift vector algebraically enter in the action. 
Hence, variation with respect to cf and B always lead to constraint equations. So, we can 
use (2.19) and (2.20) constraint equations to eliminate 4>i and Bi from the action and use 
background equations (2.16) and (2.17) to obtain a second order single variable action in 
terms of ipi. Further, writing the action in terms of Mukhanov-Sasaki variable ‘u’, 

(2)^7VC ^ 1 y rf4^|^/2 _ ^2 ^ij ^ y^2| (2.31) 

where v,z and Cs are dehned in equation (2.23). We can vary the action (2.31) with respect 
to V to obtain equation of motion of v, which is identical to the equation (2.22). Hence at 
first order, order-by-order Einstein’s equation approach and reduced action approach lead to 
identical result. 

Similar procedure may be followed to obtain reduced second order single variable equa¬ 
tion of motion. Fourth order action does not contain terms that have time derivatives of 
4>2 and B 2 . Hence, as in the first order, one should be able to substitute (f >2 and B 2 in the 
fourth order action to obtain a reduced action in terms of (/? 2 - Malik et al showed, for canon¬ 
ical scalar field under slow roll approximation, that the single variable equation from both 
approaches are same[54]. Similarly, since in the case of non-canonical scalar field, equations 
of Lapse function (/)2 and Shift vector diB 2 are (1) identical for both the approaches and (2) 


- 8 - 


are constraint equations, reduced single variable form of the equation of motion should also 
be identical. In Appendix B, we give the reduced single variable action as well as equation 
of motion in terms of V 2 ’ for non-canonical scalar field in Power law limit. Hence, both 
approaches give the identical results up to second order. 

At the first order, equations of motion are linear in first order variables. In higher 
order, only the highest order perturbed variables appear linearly, where the lowest order 
perturbed variables contribute non-linearly to equations of motion. For example, in second 
order, equations are linear in second order variables ip 2 , 4>2 and B 2 but are quadratic in first 
order variables and Bi. Hence, as pointed out in [50], it does appear that obtaining 

equations of field variables 0, B and at higher order from the two approaches may not be 
identical and thus the reduced form of the single variable equations of motion of the two 
different approaches may differ. However, instead of the non-linear form of the perturbed 
action, reduced single variable equations of motion at second order obtained from both the 
approaches are identical, hence we can generalize that at every order, in the case of non- 
canonical scalar field, both approaches lead to identical result. This leads to the following 
question: 

Why Appignani et al[50] obtained different equations of motion from two approaches? 

In the simplified model proposed in Ref. [50], authors have assumed that the homoge¬ 
neous universe is filled with matter fluctuations with no Lapse function (f and Shift vector 
diB. They have shown that in this simplified model stress tensor, energy density and pressure 
are not identical for both the approaches. Note that, since there are no metric fluctuations, 
left hand sides of the perturbed Einstein’s equations are zero. This leads to five equations 
(Tg = 0, = 0 and equation of motion of scalar field) for a single perturbed variable 6ip, 

which lead to the inconsistency of the simplified model, hence the ambiguities. Another way 
of looking into this is the following: in the action approach, one can obtain the Hamiltonian 
(momentum) constraint of the system by varying the action with respect to 4> (B). Since in 
this simplified model, both are not present, this leads to inconsistent results. 

This leads to another important question which we address in the rest of the paper is: 

For what theories of gravity and matter field, the two approaches lead to identical single 

variable equation of motion? 

To answer the question, let us look at the procedure of conventional gauge invariant cos¬ 
mological perturbation theory, which is based on two things, first, to obtain gauge invariant 
variables and second, to obtain a single variable action/equation of motion in terms of gauge 
invariant variables. Gauge invariant variables are model independent (if the background met¬ 
ric is unchanged), i.e., these are same for canonical, non-canonical or Galilean models so that 
we can always remove two variables out of five by using gauge conditions and define suitable 
gauge invariant variables. At each order, we start from five perturbed variables {4>, B, ip, E 
and ip). The gauge choice helps to remove two variables. Garefully choosing a gauge (in our 
case, E = 0 and V’ = 0) at any order could fix the gauge issue and reduced variables will 
coincide with gauge invariant variables [53]. So all equations in terms of those variables also 
become gauge invariant. 

Obtaining a single variable action/equation of motion depends solely on gauge fixing 
(the procedure discussed in the above paragraph) and two constraint equations which differ 
from model to model. If Lapse function N {(p in perturbed case) and Shift vector Ni {diB in 
perturbed case) remain constraints for any models, i.e., those functions algebraically enter 
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into the action then equations of motion of Lapse function and Shift vector contain no time 
derivatives of them, we can always eliminate them from action/equation of motion to get 
a single variable action/equation. This helps to reduce the degrees of freedom to one and 
we can write the action/equation of motion in a single variable form. However, if (j)i or Bi 
or both become dynamical i.e., if the action contains terms containing time derivatives of 
Lapse function and/or Shift vector such that equations contain double time derivatives of 
those variables then it is not possible to substitute those variables in the action or in the 
equation of motion of the scalar held and the method fails. We refer the constrained nature 
of Lapse function and Shift vector as ‘Constraint consistency’ condition. If it is satished then 
the whole method of gauge invariant cosmological perturbation theory will work. In the next 
section, we test the ‘constraint consistency’ condition for several models that are used in the 
literature. 

In fact, the whole exercise may be done in terms of Lapse N, Shift Ni and scalar held ip 
without applying any perturbation theory. From Hamiltonian theory of General relativity or 
from Einstein’s equation we obtain constraint equations, i.e., Hamiltonian and Momentum 
constraints which are functions of {N, Ni,'yij,p,p',dip) in which Lapse function and Shift 
vector are constraints. If out of four constraint equations (1 Hamiltonian equation or 0-0 
Einstein’s equation and 3 Momentum constraint equation or 0-i Einstein’s equation), we can 
solve and extract four quantity, one Lapse function and three component Shift vector then we 
can substitute those back in the action or in equation of motion of the scalar held. Unfortu¬ 
nately, GR equations are so highly non-linear that analytically solving constraint equations 
for Lapse function and Shift vector and obtaining a single variable action or equation is very 
difficult. Perturbation theory helps to simplify those equations so that we can invert those 
equations in terms of Lapse function and Shift vector and obtain a simplihed solutions of 
constraint functions. 

3 Specific models 

In this section, hrst we start with well known models of inhation within the framework of 
general relativity and then move to modihed gravity models. To check for constraint con¬ 
sistency we follow action formulation, write down second order action in terms of perturbed 
variables and identify terms that contain time derivatives of Lapse function and/or Shift 
vector. 

3.1 Minimally coupled Galilean field 

We start with the model with derivatives of metric in the action, 

S^ = I (fx^g G{X, p)np (3.1) 

which has been proposed by Kobayashi et al[44, 45] where ‘D’ is defined by equation (2.2). 
Since ‘D’ contains time derivatives of metric as well as matter, it is not obvious whether the 
action can be expressed in a single variable form. After partial integration, the second order 
matter action becomes as follows: 
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(2) cG _ 

'~'m. 


J d'*a;|5Gx0i<^i<^o^ + —GxV’o'/’iV’o ^Gxo^'(piVo^^ ~ x<P\'fi'iVo — x<t>i‘Po‘Pi‘Po + 

9Gx4>i‘fW‘fo^a.~^ - GxS^^diBidojBiip'o^ - ^Gx<5‘''a«Bi9iBi¥’oVo^ + ^GxS^^a'diBidjBup'o^a-^- 

3Gx</>i/iVo" + GxV’Wi'' + 2 Gx¥’o¥>i/i' - 3Gx¥^oaVi''a“^ - 2GxS^^^'od^Bidj^i^'^- 

GxS^^diy^idojBup'o'^ - 2GxS^^d,Biaoj‘Piv'o^+3GxS^^a'a,Bidjipi^'o^a-^ - 2GxS^^v^'oai^iaojipi + 

Gx^^^‘Po<^'ai‘piajipia~^ + Gx^'^aiBiaj(i>i‘fo^ + Gx^^^ airpiajipiipg^ — 

5Gxx‘Po<t'Wo^a<--'^'l + 5Gxx«'</>?V3o®«^“^^ + 7Gxx4>i‘PWdv'o^- 8Gxx<^i‘/’iaVo'*“^“^' + 
Gxx<Pi‘fi‘p'o‘^o,^~^^ + -Gxx'5*^9i-Bi9jSi93QVo"^a^“^^ — -Gxx^^^ o-' diBiajBupQ^ a^~^'> + 

- ^GxxV’oVo^V’i^a^”^’ + ^GxxaVo^V’i^a^"®' - Gxx'Piv^iVo^^^”^' + Gxx'S'^SiBi^jV’iV’oVo^a^'^'- 
Gxx^^-’a'SiSiajVJii^o'^a^"®^ + iGxx<5''’9i‘Pi9j¥’i¥’oVo^“^“^' - ^Gxx<5‘^a'9i¥’i9i¥’i¥5o^“^“®'- 
Gx'K'^lVo^Vl - SGxy'^iV’oVo^V’i + SGxy^iaVo^a^^V’i + 2Gxy‘(2o‘^iV5oVi - 3GxyV’iaVo^“~^¥5i + 
GxY'f'l'f'o^^i + -Gxxx‘p'd‘PiV()^<^^~‘^^ - -Gxxxo,'<l>\'f'o^— Gxxx4>ivWiv'o^+ 
Gxxx</>iV5'iaVo®<i^~®^ + -GxxxV’oVo'^'/’i^a*'”'*^ — -GxxxaVo*¥’i^“*'~®^ + -GxYY^d'f'o^vi — 
-GxYY<^''f'o^V\0,~^ + Gxxy0i¥’oVo"^“*'~^Vi — Gxxy</>iaVo^“^~^Vi - Gxxy¥5iV’oVo^“^~^Vi + 
Gxxy¥’'i«Vo^“^“^Vi + _ ^GYtpiv'o^a? + GY(t>i‘pWi<^^ “ ]^Gy^'i^cl^+ 

GyS'^ ip'o^iBidjipia'^ + iGy(5*29i(^i9j^ia^ + j^Gxy<^i‘/?o"^ “ ^Gxy^iVj'ii/^o^ - iGxy<5*'’9iBi9jBi¥5o"^+ 


-Gxy'Po^v'i^ — -GxYS'‘^aiBiaj(piifiQ'^ - -Gxy^'^^ ai^pidjtfup'g"' + -Gyy( l>ivd^o,‘‘ipi - GYY>pd‘Pi<^‘‘vi 

- -GxxY<Pi‘Pd^cL^~^^ + -Gxxy 0 i‘^i¥’o*“*'~^^ — -GxxYvd'^‘p'i^o,^~^^ - -Gyyy'p'o^ViO-^- 
-GxYY4‘i¥>d^‘fi + -G xyy’pWo^ ‘Pio.^'^ 


„'3 


2 


v 2„2. 


(3.2) 


where Gx = 9xG, Gy = d^G for background and so on. The derivatives of the constraints 
(i;i)i) appear linearly in the above reduced action (those terms are highlighted in the above 
expression). However, by performing partial integration one can rewrite these terms as terms 
proportional to (j)i e.g., first terms in the action can be written as —|9o{Gx</^o^}0i- So, 
although the action contains time derivative of Lapse function but it is reducible to action with 
no time derivative of Lapse or Shift, hence, variation of these terms do not lead constraint 
inconsistencies. 

3.2 {Lp.\(p'^{\Lnp}‘^) model 

Let us consider the following model where the matter action is given by 

Sm = j d'^Xy/^ip.xV^'^iOipV (3-3) 

and is minimally coupled to gravity^. Expanding the matter action to second order, we get, 

^Note that, our motivation is only about the consistency of the method discussed in the first section for 
different models, not to check its physical observational viability or any other problems such as Higher order 
Ostrogradsky’s ghost 
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(2)c _ 

— 


j d^x{-^4>1(fiQ^ipQ^a - YOa'tpQtplipQ^a ^ - 144ii(f>[ipQ(fiQ^a - lOtpiS^^tpQdijBiipQ^a ^- 

70<Pl(pQ'^a'^a~'^ — 28ij)i<j>ia'a~^ — 20<piS’^^a'dijBup'g'^a~^ + 10(/)i</JgViVo^“~^ ~ ‘Podjupi^p'o^a~^ 

+ 204 >ia'ip"(pQ^a~^ + 60<j>iip'ia' ipQipQ^a~^ — 124ii5'‘^ a' djitpiip^^ a~^ + 80(/>iv5j</Jg®a'^a“'* + lO0i(poi^i</5o 
+ 6(/)'iip'i(poVo^“~^ + GS’’^ V’i‘PodijBnpQ‘^a~'^ + 16<f>[ip'ia'ipQ^a~^ + 166’’^ tp[a'dijBupQ^ a~^ — 4:ipQip[tf>Q(fi'{a~'^ 
+ 45^^ tfQipitpQdjHpia~^ — 12tpia'a~^ — 12<fiQa' (fig a~^ + 12S‘^ <p'ia'dji<pi(fiQ^a~^ — 24 <Pq^ <p'i^ a'^ a~‘^ 

+ ^S^^diBidjBi^'g^ip'^^a-^ + 105^^a'diBidjBlip'll+ lOS^^diBidjBip'g'^a'^a-'^ - ipi Vo^a"^+ 

25^^ Pgdj Bidipipg^ a~^ + + 86'^^ dj Bidipipg^ a'^ a~^ + dip\djpipg‘^a~‘^ + 

4S’^^ p'ga' dipidjPipQa~^ + 4S'^^ dipidjpipg^ a'^ a~* + 46’'^ diBipgdjopip'o^ a~'^ + 2(l>[p'{pQ^a~'^+ 

45''^ a' diBidjpipgp'g^ a~^ + 25'‘^ dipipgdoj Bipg^ a~^ + 25'‘^ p”dijBipQ^a~^ + 8S''^ a' diBidjgpip'g^ a~^ — 
25'^^ ip'idjiPip'o^ a~'^ — 25''^ dij B\dikp\p'g^ a~^ + 2S‘‘^ dicpidjpipgPg’^ a~^ + 86"^^ diBidjpipQ^a''^a~^+ 

4S^^ a' dipidoj Bipg^ a~^ + 4<5*'’a'9i0i9jipi(pQ®a“^ + 25'^ diB\pgdgjB\p'g^ a~^ — 25^^ diBidjcpipgPg^ a~^ — 


' ■ 1^-2 


4>'^ ipg a ^ + 4&''^a'diB\dgjB\ p'q a ^ — 45’’^a'diBidjCpip'g a ^ — 28'^^ij}\dijB\ p'q 
5^^5’-^dijBT,dMBip'g^a-^ - p'g^pfa-^ + 28^^p’ldi.pip'g^a-^ - dup^dikPip'g^a-^} 


( 3 . 4 ) 


Following points are worth-noting regarding the above action: (i) Unlike Galilean scalar 
fields, the above action contains square of the derivative of constraints (0^^) [the term is 
highlighted in the above expression], (ii) In case of Galilean, it was possible to rewrite it 
as boundary term, in this case it is not possible. This implies that using reduced action 
approach, we cannot obtain a single dynamical equation. Hence as discussed in section (2), 
the constraint consistency is not satisfied. 

3.3 model 

As like the previous subsection, the following action also contains higher time derivatives of 
constraints. Expanding the matter action 

Sm = j d'^x (3.5) 


to second order, we get. 


= J — ^‘P^V’o^‘P'o^ o, ^ + 3ba'(pgipiVo^o, ® ^ ^-|-14())i</)ja'ipo^<i ^ — 

10<pi5'^^a'dijBiipg'^a~^ + 10(/)iipgViVo^“~^ ~ lOcpia'p”p'g^a~^ — 30<j>ip'ia'p'^pg^a~^ — 6<t>i5’'^a'djipipg^a~^ 
+ 80(j>ip'iPg^ a'^ a~'^ + 104 >ipgPiPg'^ a~'^ + G(f>ip'iPgPg'^a~^ — 8<j>'ip'ia' p'g^ a~^ + 85^^ p'^a' dijB\p'g^ a~^ — 
4p'gp'^pgp'{a~'^ + &p\a' p'Ip'i^ a~^ 4- Qp'ga'pgp'-^a~^ + Q5''^p'^a'djipip'g^a~^ — 24p'g^p\^a'‘^a~'^+ 
^&*^diBidjBip'g^p'i^a-'^ - bS^^a'diBidjBxp'^p'g^a-^ + 105^^diBidjBxp'g'^a'^a-'^ - </?; Vo^“"^ + 

26'^^ Pgdj Bidipipg^ a~^ — 48^^^ a' dj Bidipipgpg^ a~^ + 88'^^ djB\dip\p'g^a''^a~'^ + 8'‘^ dipidj pipg’^ a~’^ — 

28'^^ p'ga' dipidjpipga~^ + 68^^ dipidjpipg^ a'^ a~^ + 48''^ diB-ipgdjgp\p'g^ a~^ + 2(p'ip'{pg^a~^ — 

28'^ a' diB\djp\pgp'g^ a~^ + 28''^ dip\pgdgjB\p'g^ a~^ — 48'‘^ a' diB\djgp\p'g^ a~^ — 48'^^di<l)\djgp\p'g^a~'^ — 
28'^ 8^^ diisB\dijp\p'g^ a~^ + 28"^^ di4>idjp\pgp'g^ a~^ + 88’’^ diB\djP\p'g^a'‘^a~'^ — 28''^a'dip\dgjB\p'g^a~^+ 
28'^ a' di(l>\djp\p'g^ a~^ + 28^^ diBipgdgj Bipg^ a~^ — 28"^^ diB\dj<l>\pgp'g^ a~‘^ — cj}'^‘^pg'^a~'^ — 
28^^a'diBidgjBip'g‘^a-^ + 28^^a'diBidj4>ip'g'^a-^ + 28*^di4>idj4>ip'o'^a-^ - 8^H'^^diuBidjiBxp'g'^a-^- 


p'g p'{ a + 28'^^digpidjgp-ip'g a ^ - 8^^8’^‘dkipidijpipg a ^ - 48'‘^a'dipidjgpipg a 


( 3 . 6 ) 
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It is important to note that the second order action contains which cannot be 
absorbed as a boundary term (the highlighted term in the above expression). Hence the 
constraint condition is not satisfied leading to the fact that the reduced action does not lead 
to single variable equation of motion. 

From the above analysis, we can generalize and apply this to any higher derivative 
scalar theory models like and that obtaining a single 

variable equation of motion or action is not possible for these kind of models. 

3.4 f{R) model 

Until now, we have considered different forms of scalar field action without modifying gravity. 
In this subsection, we consider the simplest modification i.e. R + aR^ while we consider the 
matter to be canonical scalar field. Since R and matter part of the action does not have any 
inconsistency, we expand R^ up to second order in perturbed variables to get, 

36S^^ + 12S^^ S’^'-a'dijBidokiBia-^ + 725*^4>{dijBia'^a-^ + 

2S3<j>i&^^a'a"dijBia-^ + 125^^ S'^‘a'dijBidki<t>ia-^ + 126^^ S'^^a'duBidoijBia-^ + 4F'»,5“aoo-Bi9ofc!-Bia"'‘+ 
245B<f,[a'aoijBia-^ + 960iF'»a"aooBl““® + 46^^ S^''dkl<Pldoij Bia-* + 36ct>'^^ a'^ a~^+ 

432ipi(f>ia' a" a~^ + 246'^^ (p'^a'dij + 432<pfa"^ a~^ + 96(piS^^ a"dij(f>ia~^+ ( 3 - 7 ) 

126^^5^’^ a'dkiBid^jcPia-^ + 45Bs'^‘dij<PidokiBia-^ + 4SB5'^‘dijcPidki<pia-^ + 24S^^cl>[a''dijBia-^- 
72S^^a'diBia''dojBia-^ + 12S^^S'‘''a"aijBiakiBia-^ + 72S*^a'aiBiajtpia"a-^ - 72(5*J9iBio'jBia"^a-® + 
245Bai4>iaj4>ia''a-^ - 125^S'^^a"aikB^a^iBia-^ 

Following points are interesting to note from the above expression; (i) R^ term contains 
time derivative of (j), that cannot be absorbed as a boundary term. This implies that the 
above action cannot lead to the constraint equation, (ii) By doing a conformal transformation 
—)■ = fl'^gfj^i,, the term containing the time derivative of (f) can be absorbed as a matter 

field and hence the constraint equation recovered, (iii) The constraint consistency allows us to 
identify that the f{R) gravity models, without conformal transformation lead to inconsistent 
dynamics. 

3.5 [ip-\(p'^{{nipY - model 

As we have shown above, certain higher derivative models do not satisfy ‘constraint con¬ 
sistency’. We have also shown which terms in the second order action spoil the ‘constraint 
consistency’. However, it is interesting to note the terms that contain time derivative of 
Lapse functions are identical. Let us consider the following scalar field action 

S = j ] (3.8) 

The second order action is given by. 
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( 2 ) 5 ^ = J d'‘a;{- 105 aVo</'i¥’o^a“^ - lOrpiS’^ifodijBiifo^a-^ - 42 <^l(/)'iaVo'‘““^ - a'dijBiip'o'^a 

G(j>iS^^ ipQdjHpi(pQ^a~^ + 30(j>ia' ip"ipQ^ a~^ + 90<j>np'ia' (fQip'o^ a~^ — 6<j>i5'^^ a'djUpiipQ^ a~^ + 

66 '^^dijBiip'g^a~^ + 24 cf>[(fi[a'(fiQ^a~^ + ipga'dijBiipg^a~^ + 4 S^^if!gip\^f:gdjiip\a~'^ — 

\8^\o! a~'^ — 18</jQaVoVi^“~^ + ip\a! djnp\ip'i^ a~^ + 155^^a'diBidjBiipg!pQ^a~^ + 

125 '’^ a' dj Bidiifiiipgifg^ a~^ + 65 ’’^ ip'ga' diipidjipnpQa~^ — 26 '^^ dupidjipiipg'^ + 66 ^^ a'diBidjtpiipgifg^a 

2 S^^ip”dijBiipg^a-^ + 125 ^^a'diBidjoV^iVo^a.-^ - 2 S'^ 4 >{dji^piipo^a-^ - 2 S'^ 5 ^‘dijBidik‘PiVo^a.-'^+ 
a' dupidoj Bupq^ a~^ + 2 S''^ a' di<f>idjipiip'g^ a~^ + 65 ^^ a' di Bidgj Bupg'^ a~^ — 66 ’’^ a'diBidjipi‘Po^o,~^~ 
26^^<t>[dijBiip'g^a-^ - S'^5^‘di,BidkiBi^'o*a-^ + 2S'^^'{dji^iip'g^a-^ - + 

4 S^^di<PidjoV^iVo'^a-^ + 2 S'^ 5 ^‘dikBiaijViip'o^a-^ - 2 S‘iai^iaj^Kp'g'^a-^ + 5 *^ 5 ^‘aikBidjiBiip'g‘^a-^- 
26 ’-^ aio‘fiajQ(pi(fig’^ a 


J 2„-2 ,. s'^5>‘^aki<fiiaij<fii:fi'o^a-^ + AS^^a'ai^iajQifiw'o^a-^} 


2„-3 


(3.9) 


It is interesting to note the following points: (i) the action does not contain terms 
having time derivative of Lapse function/Shift vector. Hence the resultant equation leads 
to constraint equation. Similarly, — 3n<y9 + 2(p’.^(p’^ip’^)] model do 

not have dynamical Lapse/Shift, (ii) From the above, one may be tempted to relate the 
constraint consistency with Ostrogradsky’s instabilities. To go about checking this, let us 
look at the zeroth order action for the matter field, i.e.. 


= J d'^xj-eaVoVo^a^ 

which, after integration by-parts, can be re-written as. 


(3.10) 


(0) C / j4 " -3 /4 y ,2 -4 /4.| 

^ ’Sm = d x{-a a - -a a ] 


(3.11) 


Hence the equations of motion will contain a'" and (I)q . This immediately signals Ostrograd¬ 
sky’s instability. 

This leads us to the important conclusion: While all ‘constraint’ inconsistent models 
have higher order Ostrogradsky’s instabilities but the reverse is not true. One can have 
models with constraint Lapse function and Shift vector, though it may have Ostrogradsky’s 
instabilities. 


3.6 Constraint consistent models withont instabilities 

In the previous subsection, we showed that identifying the terms in (3.4) and (3.6) that 
contains higher derivatives of Lapse function, we can remove these terms by combination of 
these terms two terms. However, we noticed that such actions suffer from Ostrogradsky’s 
instabilities. The term that leads to Ostrogradsky’s instability is at zeroth order 

action (3.11). In order to cancel such term, without involving any higher derivatives of f, 
one needs to introduce non-minimal coupling of the kinetic term of the scalar field with Ricci 
scalar. It is easy to show check that the term — |(y?;cKy9’“)^i? exactly cancels this and the 
resulting background action becomes, 




(3.12) 
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Similarly, [(/j.({□y?}^ — SDf/? ip'-^ip-aip-^u)] does not have 

any constraint inconsistencies as well as instabilities. 

The following points are worth noting regarding this: 

i. We have arrived at the action — ^{ip-^a(p'°^YR] by using the 

condition that the action does not have time derivatives of Lapse function and Shift 
vector and later using the additional condition that Ostrogradsky’s instability does not 
arise. 

ii. In a different manner, Nicolis et al[46] and Deffayet et al[47] have come up with the same 
action, only with condition of removing Ostrogradsky’s instability. Here we have veri¬ 
fied that those models will result in consistent dynamics with ‘constraint consistency’ 
and lead to single variable action as well as equation of motion. 

iii. Similarly, we find that the third order derivative model is ip-^yip’^’^-\- 

2(p]^(p':^(p'^—QG^^(p'^(p-a^-p)] which, again has been derived in [46] and [47]. This model 
is also constraint consistent and free of Ostrogradsky’s instability. 

4 Conclusions 

In this work, we have revisited the two approaches in cosmological perturbation theory — 
order-by-order approach of the Einstein’s equation and reduced action formalism. Equiva¬ 
lence of the two approaches were not clear since Gravity equations are highly non-linear. In 
this work, we have established that equations arising from order-by-order approach of the 
Einstein’s equations and those from the action formalism are equivalent for canonical as well 
as non-canonical scalar fields up to second order. These results may easily be extended to 
any model in Gravity models at any order. 

To compare both the approaches. We have identified a ‘Constraint consistency condi¬ 
tion’ where the constrained nature of Lapse function and Shift vector are studied. We have 
shown that, in order to obtain a reduced equation for both of the approaches, ‘Constraint 
consistency’ relation has to be satisfied, i.e., those variables should appear in the action al¬ 
gebraically, and no non-reducible partial time derivatives of Lapse function and Shift vector 
should be present. In other words, equations of motion of Lapse function and Shift vector 
should not contain second order partial time derivatives. We then investigated the higher 
order derivative theories of gravity and found that models which satisfy the constraint condi¬ 
tions can be applied to the conventional perturbation theory where we express all equations 
in a simplified form with a single variable (Curvature perturbation TZ or Mukhanov-Sasaki 
variable) equation of motion. One common problem with higher order derivative theories 
is that, they may have Higher derivative equations of motion which can give rise to Os¬ 
trogradsky’s instabilities. We showed that all the models which do not satisfy ‘Constraint 
consistency’ conditions suffer form Ostrogradsky’s instabilities. However, we also find that, 
there exist some models which satisfy constraint conditions though they have instabilities, 
i.e., the action can be reduced in a single variable form but the single variable equation of 
motion will contain higher order time derivatives. The method we have proposed here is fast 
and efficient and would be useful for inflationary model building. 

We have also constructed some higher derivative models in such a way that those mod¬ 
els should satisfy constraint consistency condition and can be free from higher order time 
derivatives. Those models have already been derived in the literature in a different manner 
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where the approach of constructing models are different. This ensures that, conventional 
perturbation theory can be used in higher derivative models which are free from Ostrograd- 
sky’s instabilities to obtain a gauge invariant single variable equation of motion using any 
approaches. We summarize the main results, 

1. The two approaches are completely equivalent. 

2. Not all models with Lagrangian density £ = f + F()f, dip, or £ = 

f y/^F{R, dp, d'^p) can be reduced in a single variable form. 

3. To obtain a single variable form, ‘Constraint consistency’ condition has to be satisfied. 

4. Models with constraint inconsistency show Ostrogradsky’s instability but the reverse 
is not true. 

The analysis here may have implications for models of quantum cosmology with scalar 
fields[55]. The quantum corrections to the matter and gravity can be modelled as effective 
stress-tensor [56]. The effective classical equation must also satisfy the constraint consistency. 
This might help to constraint the quantum cosmology models [56]. 
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A Coefficients of second order equation of motion of non-canonical scalar 
field 


Cx =P2 - - 4'2 Po - - 4>2p'o + 2'Hp2 + ipodi - -I- 6(j)i4>Wo- 

mdip'i - 2H(I)2 Po - PoV^R2 - 2p[V^Bi - AS^^diBidjp[ - di(t>id^pi - 25^^d^pidjB[ + 
GHp'odl + 2(I3 ip'qS/^Bi + - 26^^ p'^d^BidjB^ - d^Bidjpi- 

S^^d^BidjBip^' - 2R5^^p'^d^BidjBi 

(A.l) 


Cxx =a - 3a - 12a (l)i4>'iPo'^ + 8a ‘^(t)iPiPo‘‘+ 

-I- ~ a~'^4>2'Hp'o + 2a~‘^p'^p^V'^pi-\- 


,-2,JI,J 2 


Aa~‘^ 5''^ p^dipidjPi + 2a~^ p^V"^ Bip^^ — 2a~^S’‘^ diBidj(l)iPo^+ 

2a~'^5''^diBidjB'ip'^ P diBidjp'ip'^ — 2a~^6^^dididjPiPQ^ + dipidjpipQ + 

2a~'^5^^dipidjB'lp'o^ + a~^'Hp2Po^ + 3a~^'H(piPo^ ~ 2a~‘^'H4>ip'iPQ^+ 
6a-^S^^p'od,BidjPip” + 4.a-‘^5^^dM^Bip'^p'^ - d^B^d^B^p'^ - 

2a~‘^'H5^^ diBidjpip'o^ 

(A.2) 
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Cxxx =‘^a - 2a + a '^'H4>2<p'o'' - a '‘(/)2V5oVo^ - 2a '‘^iv3it/?o‘*+ 


_4 / // / u —4o / / / I —4 / // / 'J O —4o / / 'J r —4o / / '^ / ■* i 

2a ipiifiipQ —a H(p 2 ‘Po +a ‘P 2 Vo‘fo — oa ncpiipQ — oa H(Pq (pi + 


11 a + 6 a "‘(/SgVo^'/ji^ + 12 a ‘^'H4 >\(p'i(p'q* - 16a "‘(/)i(/?'i(^oVo'^+ 




—4n 


.-4. 


,JfJ 3 


(A.3) 


+ 2a-^'HS^^ diBidjifiip'^ + a-^nS^^di^pidjifiip'o- 

a~'^6'^^diBidjBiifQip'f^'^ - 2a~'^5'-^d^BidjipiipQipQ^ - a~'^6'^^diipidjipiipQip^^ 

Cxxxx =a~^'H(j)\ip'Q + a~^'Hip'o^- a“®(/?o^iV5o^ - ar^^'W^v'\- 

2a“®'H^iV3i(Po® + 2 a“®(/)i(/?'i(^oVo^ 

Cxxxip =a~^<^\ip'Q + - 2ar^(i>xp\p'^ + 2a~‘^'H(t>ip'o^pi - 2o“‘‘(/)i(/?oVo^'/5i- 

2a~'^'Hp[pQ'^pi + 2a~'^p[pQPQ^pi 

CxXip =a~'^(l)2p'o^ - a~^p2Po^ - 5a~'^(j)lpQ^ - 4a“^v5o^'/5i^ + 8a~‘^(j)ip[pQ^+ 

2a~‘^(j)[pQ^Pi - or‘^p'^p'QP2 - ‘^a~'^p'[p'^px + 8a“^()ii(/?oVo^'Pi + 

a~^ 5^^ diBidjBip^^ + 2a~^6^^ diBidjpipQ^ + a~‘^6^^ dipidjpip'^ — 6a~'^ p'qp'iPqPi + 
a~‘^'Hp'^^ P2 — 2a~‘^H(j)iPQ^ Pi + 2a~'^'Hp'ip'Q^ pi 


(A.4) 

(A.5) 

(A.6) 


/-i o —2i n —2 / /3 —2 // / ^ 2 I —2o/ /2 

Cxx<p<p=2a (piPoPi-2a p^pg pi - a pgPg p^ + a Bpg Pi 


2,„2 


3,„2 


(A.7) 


Cxv =‘/3'l^ + ‘Pot/32 + V3 oV2 + 2p”pi + <j)lp'g^ - 2(l)ip'gp\ - 2(l)lPgPl- 

S^^dipidjpi - 25'‘^dijpipi - 2(j>'-^p'gpi + 2'HpgP2 + 43-Lp\pi- (A.8) 

^^.(pip'opi - 25''^p'gdiBidjpi - 25'^^p'gdijBipi 


Cxipifi — ‘PoV’i A T’o ^^2 + 2pqPiPi + 2'HpgPi (^-9) 


/ 2 2 

^Xipipip '^1 

(A.IO) 

= a^(j)2 — a^cjil + a^6^^ diBidjBi 

(A.ll) 

Cipip = P2 + 2(X^ (j)\p\ 

(A.12) 

n — ^2_2 

^ppp — ^ Y^l 

(A.13) 
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B Second order single variable equation of motion for non-canonical scalar 
field for Power-law infiation 

Fourth order action for non-canonical scalar field, after partial integration is, 


( 4)5 = J d*x(3<f>i(l}2S"^a'dijBiK-'^a+ ^Px+ ^P<t>W + ^Px<t>2S^^diBidjBi^p'o'^a^ + ^P(l>25^^diBidjBia^- 

-Px<t’2<l>iPo^a^ - 3P<p2<j>iO-‘'‘ — -4>2S^^dikBidjiBiK~^+ -ip2S^^S^‘dijBidkiBiK~^a^ - -Px4’i<t>2PoV^iO?+ 

-Px<l>2VoP2°‘'^ + -PxS'^WW2diBidjBia^ — -Pxip'o<fi24’iO-^ + -Px4’2‘fi'i^+ -Px4>ivW2°‘'^ - -Px9'2^o? — 

4 4 4 4 2 o 

]^Px<t>2&''^'p'odiBxdj^ra^ - ^Px<l>iS^^ipodiBidjip2a^ + ]^Px&''^p'idiBxdj^2a^ + 


9 , 


-PxS^'^ diip2djip2a^ — -PxX<p24'l‘P'o* + 3PxX<i>\<f>2V’Wo^ + -PxXP2<l>i‘fio'^ — '^Px X + 


/ J.2.J 3 


-Pxx<l>iS^^diBidjip2‘Po^ + -Pxx<l>iS^^diipidjtp2‘Po^ + -Pxx<f>2‘Po‘^ + -Pxx4’2S^^diBidjBiifig'^ - -Pxx<t>2‘P2Vo^ — 

2. 2 o 4 4 

Pxx<t>2Vo^‘Pi'^ + -Pxx‘p2S'^^ diBidjipiipQ^ + -Pxx<t>2S^^ diipidjipiip'o^ — -Pxx^^^‘f29iPT-djBnpQ^+ 

-Pxxv'oVWl’^ — -PxxS''WldiPldj‘P2Vo^ — -Pxx&^^ V^0‘p'ldnpidj(fi2 + -PxxVo^‘P2^ - -Pxx^^^ V22diPldjipiipQ^ 

4 2 2 o 2 

- - Pxx^^^ <PQip2dnpidj(pi + -PYY<P2a* + -Pxy4>2'p'o^O?V 2 + -PxYS^^diBidjBYPQ^a^tp2 - -PxY<Pl'p'o^O?V2 

4 o 4 4 2 

- PxY<l>l(t^ 2 Po^a^<fil + -Pxy4>WoP'i°-'^‘P 2. + -PxY<t>2<fio‘fil<^'^Pl + -PxY<l>l‘Po'fi2°'^‘fil - -PxY‘fioP 2 <^'^‘fi 2 — 
-PxyPi^o-^'P 2 — -PxY'fii<fi2a'^ipi + -PxY^^^ PodiBidjipia^(f2 + -PxY^'’^ v'QdiBidjip2a^(pi + ■^Pxy^''^ diipidjipia^(p2 + 

-PxY^''^dnp\dj‘^20?V’l + -PxXx4>2<t>lPo^O'^~^^ - -PxXX‘P24>1‘Po^O,^~^^ — -PxXx4‘l4'2V'lP'o^a!'~^^P 
-PxXX<t>l‘f'l‘fi2Po‘^<^^~^^ + -^PxXX<t>2‘f'o*‘‘p'l^a^~^'’ — -Pxxx<fi2‘fi0^'fil'^<^^~^^ + -PYYY‘fiia^‘P2 + -PxXY<P'l‘Po*‘fi2 + 
-PxXY<t>l4>2‘Po*Pl — -PxXY4’lPl‘P'o^‘f2 — -PxXy4>1V2Vo^‘P1 “ -PxXY4>2‘PlVo^Vi + -PxXY‘fo^Pl^<fi 2 + 

-PxXY‘fi'i‘P 2 ‘Po^Pl + -PxYy4>1‘Po'^ 0?‘P1‘P2 + -PxYY<l>2V'()^‘Pl<^^ — -PxYYV’oV'iO-'^ PlV 2 — -^PxYY‘foP 2 ‘PlO-^ + 

-Px<t>i^’'^ diifidjifi2a^ + -PYY<l>io-‘^‘fiiP2 + -^Px<I> 25’^^ dnpidj(pia^ + -^Py 4’20 .*‘<^2 + -^PYY4>2‘fil<^^+ 


-PyS^^ diBidjBia*ip2 — -Py4‘\o-*P2 — -JV</'i 02 «^V’i) 


(B.l) 


We can eliminate (j) 2 , (pi and Bi using constraint equations. While this is possible for 
canonical scalar field, it is highly non-trivial for non-canonical scalar fields. Hence we consider 
Power law inflation to reduce the action in a single variable form. 

For Power-law, we have 


a = ao{-rif 

(B.2) 

v-^' - ^ 

i-ri) 

(B.3) 

a" _ /l(l-/3) 

a {—vY 

(B.4) 

<^o = ‘Pc{-riY 

(B.5) 

ipo' = ipcai-rj)^^ 

(B.6) 
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(B.7) 


+ P^a^) 

-2—+n‘^ = KPoa^ (B.8) 

a 

For Power-law inflation, P{X, </?) = ^ -|- X). Solving for C, pc and a using above 

equations, 


(X — 1 -j- /? 

-6/3(1-/?) 


C = 


^Oc — 


K{l+I3f 
^ do 


'1-2/3 


31 + /3V 1-/3 


Using first order Einstein’s equations. 


where. 


'^^■V 2 V 

(/)i = Fspi 


1 — /3 1 + /3 

1-2/3 ao(-r/)i+/^ 




V^Bi = 


do 


1-/3 


1-2/3 


(1-3 


V^-Bi = Pipi -|- P2P1 


1 + /3 


Pi + 


(-r/)2+^ (-^)i+/3 


7 ^ 1 } 


Fi = - 
Oo 

F 2 = - 
do 

F3 = - 
do 


' 1-/3 

1-2/3 

^ 1-/3 

1-2/3 


(1-3/3) 


(1-3/3) 


1 + /3 

(- 77 ) 2 +/^ 

1 

(— 77 )!+/^ 


'1-/3 1 + /3 

1-2/3 (-77)!+/^ 


(B.9) 

(B.IO) 

(B.ll) 


(B.12) 

(B.13) 

(B.14) 

(B.15) 

(B.16) 

(B.17) 

(B.18) 


Similarly using second order Einstein’s equations. 


9i(/)2 = CidiP2 + C2pidipi + CsSFdjpidikBi + C^d^djBidikBi -h C^p^dipi (B.19) 
</>2 = Cip2 + ^C2 pI - S^^diBidjBi + Q (B.20) 

Q = d-^^CsSFdjPidikB^ + C^p'Api} (B.21) 
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where 


Cl 


C 2 


C 3 


[3 I 1-/3 1 + 13 

V 2y 1-2/3 ao (-r/)i+^ 


3 (l + /3)(l-/3)(3-14/3+ 7/3^) 1 

2al (1 - 2/3) (_r;)2+2/3 


/3 / 1 - /3 (1 + /3) 1 

V2y 1-2/3 ao/3 (-r/)/^ 


C4 = -2 

9 (1-/3)2(1-3/3) 1 

2al /3(l-2/3) (-r/)i+2/3 


Multiplying the fourth order action (B.l) by |(1 — 2/3)^, 


(B.22) 

(B.23) 

(B.24) 

(B.25) 

(B.26) 




(Ahp 2 + A 2 d\ 2 di +2 + ^31^2+2 + + 9V2 (^Biipidiipi + B 2 +idiBi+ 

B3+'idiBi + Biiip'idi^p^ + (/92 (Di^\ + D2+1+1 + D3Q + Diiip'i + D^d''Bidiipi + 

2 V (B.27) 

D^d^ifidiipi) + ip2 {Eiifl + E2{V^Bif + E3^^^Bl^^JBl + E^+i+'i + E^+'i + E^Q 
+ Ejd^'Bidiipi + Egd'-ipidiipi 


where, 


4li = - 
A2 = — 
A3 = - 
A^ = 


3/3(1-/3)(l-2/3)(l-3/3) 

k( —T/)^ 

/3(l-/3)(l + /3)(l-2/3) 

K{—riY 

2/3(1-/3)(l + /3)(l-2/3)(l-11/3) 

K{—r])^ 

3/3(1-/3)(l + /3)^(l-2/3)(l +3/3) 

k(-?7)4 


(B.28) 

(B.29) 

(B.30) 

(B.31) 
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(B.32) 


36/3(1-/3)(l + /3)(l-2/3)(l-3/3) / 1 -/3 
\/6ao k(— y 1 ~ 2/3 
4/3(1-/3)(l + /3)(l-2/3)(l-11/3) 

12/3(1-/3)(l-2/3)(l-3/3) 
k{—viY 


^ 48/3(1-2/3)^(l-/3) I 1-/3 

^ \/6 ao k(-77)2+/5 ^1-2/3 

^ _ 9^(1 - /3)(1 - 2/3)(l + /3)(1 - /3)(3 - 5/3 + 5/3^ - 83/3^) 

\/6 ao «:(— 

^ _ 36/3V(l - /3)(1 - 2/3)(l - /3)(1 + /3)(7 - 17/3) 

\/6 ao k(—7/)^+^ 

^ -12ao/3V(l-/3)(l-2/3)(l-2/3)(l-3/3) 

\/6k(— 77 )^“^ 

^ _ 72/3V(l - /3)(1 - 2/3)(l - /3)(1 - 2/3) 

^ \/6 ao k(—77)2+/3 

_ 12/3(1-/3)(l-2/3)(l-3/3) 

^ k(— 77)2 

^ 24/3V(l-2/3)(l-/3)(l-2/3)(l-/3) 

\/6 ao k(— 77 )^+^ 

^ 9v/(l - /3)(1 - 2/3)(l - /3)(1 + /3)"(3 - 19/3 - 3/3^ - 77/3^) 

\/6 ao k(— 77 )®+/^ 

^ ^ 2ao v/(l-/3)(l-2/3)(l + /3)(l-2/3) 

\/6 ^(— 77 )^“/^ 


^ 2 ao V(l-/3)(l-2/3)(l + /3)(l-2/3) 

^ \/6 k(— 77 )!“^ 

^ ^ 72/3V(l - /3)(1 - 2/3)(l - /3)"(1 + /3)(3 + 11/3) 

\/6 ao k(— 77 )^+^ 

^ ^ 18/3v/(l - /3)(1 - 2/3)(l - /3)(1 + /3)(7 - 17/3) 

\/6 ao k(— 77 ) 3 +^ 

^ 12 gp / 3^{1 - /3)(1 - 2/3)(l - 2/3)(l + /3)(1 - 3/3) 

\/6 /i(— 77 ) 3 +^ 

_ 4/3(1-/3)(l-2/3)(l + /3)(l-11/3) 

^ k{ — ViY 

^ 18/3V(1 - /3)(1 - 2/3)(l - /3)(1 + /3)(1 - 3/3) 

\/6 ap k(— 77 )^+^ 


(B.33) 

(B.34) 

(B.35) 

(B.36) 

(B.37) 

(B.38) 

(B.39) 

(B.40) 

(B.41) 

(B.42) 

(B.43) 

(B.44) 

(B.45) 

(B.46) 

(B.47) 

(B.48) 

(B.49) 

(B.50) 
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After taking partial derivative, 


J d'^x (^Ahp2 + A2d'"ip2dnp2 + A^ip^ + Lp2 {Giip\ + E2{V‘^Bif + E^d'-^ BidijBi + 

G2 Pip'i + + Gi^pip'l + G^p'ip'l + G%pV‘^pi + G 7 (^^V^(^i+ 

+ G(^d^ p'ldipi + Gio9*(/3i9ji?i + Gud^ p'ldiBi + G 12 Q + GiaQ^ 


(B.51) 


where, 


Gi = £;i - B2F1 - D[ (B.52) 

G2 = £'4 ~ {B2E2 + B^Ei) — [D'2 + 2Z)i) (B.53) 

G3 = E5- B3E2 - {D2 + (B.54) 

G4 = -1)2 (B.55) 

G5 = -21)4 (B.56) 

Ge = -Bi (B.57) 

G7 = —B4 (B.58) 

Gg = ^8 - - (Z)^ - D3E3) (B.59) 

Gg = — B4 — 2Z)g (B.60) 

Gio = Et-B2- (Z)' - 2ZZZ)5) (B.61) 

Gii = —Z?3 — ZZg (B.62) 

Gi2 = Eq — D'^ (B.63) 

Gi3 = -Ds (B.64) 


so the equation of motion of scalar field for power law K-Inflation is, 

2Aip^^ + 2A[p^2 + 4 I 2 VV 2 + A^pI = Gipl + E2{V^Bif + E^d^^ BidijBi+ 

G2PiP\ + G^p'i + Gi^piPi 4“ G^Pip'i + Gg(^V^(/?i + G7</5^V^(/?i+ (B.65) 

G^d^pidipi + G<^d^ ip'idipi + Gio( 9 *(/Jic)jZ?i + Gud^ p'ldiBi + G12Q + GisQ^ 


oq dependency in the action as well as in the equation of motion can be resolved by 
rescaling pi —> (p 2 —> aop 2 - 
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